Flat slab deformation caused by interplate suction force

In Table S1, we provide information on the events we used for the receiver functions (RFs). The time-axis RFs are shown in Figure S1. With Figure S2 and Text S1, we derive the relationship between the change in estimated Moho depth and the change in the Vp/Vs ratio used in the reference model to convert Ps arrival time to depth. With Figure S3 and Text S2, we rationalize our observation of positive PpPp phase in the RFs. In Figure S4, we show the density of the seismicity along A-A’ profile, in which two clusters are shown, possibly related with the bending of the slab.